Temperature-dependent logarithmic corrections in the spin-1/2 Heisenberg chain. 



OS 

o> 
os 



Victor Barzykin 
National High Magnetic Field Laboratory, Florida State University, 
1800 E. Paul Dirac Dr., Tallahassee, Florida 32310 

We obtain the logarithmic corrections to the dynamic response function and NMR Ti and T2G 
rates in the spin-1/2 antiferromagnetic Heisenberg chain using perturbative renormalization group 
in the leading irrelevant operator. The result is compared with NMR experiments in Sr2CuC>3. 
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Quantum spin chains have attracted considerable in- 
terest for a long time, due to both unconventional physics 
of the ID materials and sophisticated theoretical meth- 
ods used in the analysis of the problem. In particular, 
it is well known Q that the antiferromagnetic spin-1/2 
XXZ chain is critical for J x > J z (the two exchange 
constants), and the spin correlators at T = decay as 
a power law with distance. The critical theory can be 
mapped on a theory of free bosons. For J x < J z the sys- 
tem acquires a spin gap. At the quantum critical point, 
the Heisenberg XXX chain, the spin correlation func- 
tion has logarithmic corrections to the free theory com- 
ing from the leading irrelevant operator, which becomes 
marginal. Thus, the asymptotic power law behavior of 
the correlator at long distances is modified (2|j|. The 
staggered component of the spin-spin correlator has the 
following form: 



<S*(r)S*(0)> = (-1) 



(27r) 3 / 2 r ' 



(1) 



where the non-universal coefficient D — l/(27r) 3 / 2 has 
been determined from the Bethe Anzats |^||| . 

Although this logarithmic behavior has been known for 
some time @,||, numerical [^|,0J|J|,0 and experimental 
jll],[l2| tests of it have appeared only recently. To fit the 
numerical data, various anzatses were used [pp| Jl0| , [l3|] . A 
consistent field-theoretical derivation of the logarithmic 
factor for the equal-time correlator in a finite-size chain 
at T = has been recently done in Ref. 0,0. In what 
follows we use a similar approach to obtain staggered 
dynamic spin correlator, up to two-loop order in the per- 
turbative renormalization group scheme. The imaginary 
part of the staggered dynamic spin susceptibility is ob- 
served in the inelastic neutron scattering experiments, 
and determines NMR T x . NMR T 2G is found from the 
real part of the spin susceptibility. 

Let us now show how the logarithms appear in the 
time-dependent staggered spin correlation function. For 
this purpose it is more convenient to work in real 
space and imaginary time, so that one can easily ap- 
ply bosonization and conformal invariance. One should 
then Fourier transform and analytically continue the re- 
sult to real frequences. At low temperatures and large 
distances we can use the continuum approximation. The 



theory can be written in terms of free bosons defined on 
a circle. In case of SU (2) - symmetric Heisenberg model 
it is more convenient to use the non-Abelian bosoniza- 
tion (|], which respects the symmetry. The action for 
the 5'C/(2)-symmetric matrix field g§ includes the Wess- 
Zumino term with coefficient k = 1. The Hamiltonian 
density takes the following form: 



H = H -8tt 2 /V3X3 l -J b ,, 



(2) 



where Hq is the Hamiltonian density for a free boson, 
3l,r are left and right SU(2) currents: 



A^ 



irfg+d-go-], Jj 



A^ 



Md+ggV] (3) 



The spin operators can be written in non-abelian 
bosonization notation as: 



Sj = (Jz + Jr) + const i(-iytr[gcr} 



(4) 



so that the correlation function has uniform and stag- 
gered terms, 

X (r, r) = (S Z ^> -> Xu (r, r) + (~l) r Xs(r, r), (5) 

where r is imaginary time. Both terms vary slowly on 
the scale of the lattice spacing, and correspond to dif- 
ferent Green's functions in the continuum theory. The 
staggered susceptibility, which is enhansed near q = tt 
wave vector, is observed by inelastic neutron scattering 
and NMR, 



Xs(r,r) cx (ir( CT z g)(r,r)tr( C r z g)(0)). 



(6) 



It is not difficult to determine the contribution of the 
free boson with radius R = 1/ V2n |ll| - the conformally 
invariant WZW model on a circle of length (3 = l/T in 
the imaginary time direction. Indeed, g has scaling di- 
mension 1/2, so for an infinite system one writes: 



<ir[g(z,z)<x z ]tr[g(0)<r z ] >= 



1 



zz 



(7) 



, where z = r + Here we have chosen a 

convenient normalization for the operator g. Here and 
below we use the units c = kB = h = 1. Making a 



1 



conformal transformation from the infinite plane to the 
cylinder, one easily finds: 

vrT 



< fr[g(^)er z ]ir[g(OK] >= 



(8) 



\J sin(7rTz) sin(7rTz) 

The result for real time is given by a straightforward ana- 
lytic continuation, r = it. To obtain x"(?: tL ')i one simply 
performs a Fourier transform. Integration over q in the 
limit lo -> gives NMR T x 0. These results are, of 
course, well known. 



In order to obtain the logarithm, one has to go further 
and do a perturbative expansion in the leading irrel- 
evant operator in Eq.(j|), and then collect the diverging 
terms in a Renormalization Group (RG) scheme P),p"4"l. 
To the first order in the leading irrelevant operator the 
correction can be easily calculated. The details of a sim- 
ilar calculation of this diagram in a finite-size chain at 
T = can be found in Ref. H]. The first-order correc- 
tion has the following form: 



S < tr[g(z,z)cr z }tr[g(0)cr z } >= 



vr 2 A T 



•v/3 sin[7rTz] sin^Tz] 



In 



sin[7rTz] sin[7rTz] 



const 



(9) 



Here Ao is the "bare" coupling constant for a theory 
defined with a cutoff at T = Tq. 

We can now sum the leading logarithmic contributions 
using the standard Callan-Symanzik RG equations for 
the staggered spin correlation function Xs{ r i T i T, A): 

[~d/d\nT + f3(X)d/d\ + 2 7 (A)]x s (r, rT, tT, A) = 0, 

(10) 

where /3(A) is the beta function for the coupling constant 



A in Eq.(|): 



dX 
d\nT 



"/3(A) 



(11) 



and 7(A) is the anomalous dimension. In Eq. (|T(J) the 
T-derivative acts only on the first argument of Xs] tT 
and tT are held fixed. The solution of Eq.(^0|) can be 
written as follows: 



Xs(r,r, T, A ) = exp -2 



A(T) 



7[A] 
0{K 



dX' F[X(T),rT,rT], 



(12) 



where Ao = A (To) is the "bare" coupling - a coupling 
at the energy cutoff scale To, F[A(T), rT, tT] is an arbi- 
trary function of the effective coupling constant at scale 
T, X(T). 

Since the coupling constant flows to zero as T is de- 
creased, one can use perturbative expressions for 7(A) 
and ,0(A) to determine the long-distance asymptotics for 
the staggered spin susceptibility. The universal terms in 
the perturbative expansion for the /?- function jlTj and 
the anomalous dimension are known, 

/3(A) = -(4tt/V3)A 2 - (1/2)(4V\/3) 2 A 3 (13) 
7 (A) = 1/2 - (n/V3)X. (14) 



Thus the effective coupling is given by: 



X s (r, t, T, A) 



4W 1 



v T , . jMA/T) + - ln[ln(A/T)] } + 0(1), (15) 



where 

A = const ■ ^e^^ Xo] T . 
Thus, we can rewrite the integral in Eq.(fi2| 

A(T) 7(A) 

ao m 



dx= 1 -^ + 1 -m x ^ 



T 



Ao 



(16) 



(17) 



/ Ap 
A(T) 



In general, one can expand the staggered spin suscepti- 
bility Eq.(|l2|) in powers of A(T): 



(18) 



OO 



m=0 



The coefficients, a n and the functions F m (rT,TT) can 
be determined from the perturbative expansion in the 
leading irrelevant operator. 

We can now RG-improve the perturbative result 



Eq.(^), Eq.(^). This can be done by expanding Eq.(|lg|) 
to first order in the bare coupling constant, Ao, and com- 
paring it with the perturbative calculations. We find: 



1 ^yinA + I m (lnA) , j x 

Xs (r, r, T, A) = , ■ f ^ 1 + ~TTTTTf\ HM*T*) M^Tz)] . (19) 

(27r) ci / 2 ^/sin(7rTz) sm(7rTz) \ 41n(A/T) ) 
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This expression is our final perturbative result, which 
has to be Fourier-transformed and continued analytically 
to real frequencies. For this purpose it is more convenient 



to work with temperature-dependent anomalous dimen- 
sion rj(T) instead of the form Eq.(ft9|) with the logarithms. 
To the same order we can write: 



Xs(r, t, T, A) 




In -) (sin(TrTz) sin(7rTz))^ n 



(20) 



where 



V(T) = 1 - 



21nf 



(21) 
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The analytic continuation of Eq.(ECf) is then analogous 



to that in case of a Luttinger liquid, which is well known 
We therefore only discuss quantities which 
by inelastic neutron scattering and 
NMR. Continuing Eq.(|20|) to real frequencies, we get: 



can be measured 



Imx{q,u) 



2^ T )- 2 / A 1 / A 

(>F^ sin( ^ (T)/2) V ln r + 2 ln ( ln r 



x Im < B 



i(w-g) , rj(T) rj(T) 



4ttT 



, 1 



B 



ifa + g) , V(T) 
4ttT 4 : 



(22) 



r)(T) 



where B(x,y) = T(x)T(y)/T(x + y) is the beta- 
function. An immediate consequence of the temperature- 
dependent anomalous dimension is that the correlation 
length acquires additional logarithmic temperature de- 
pendence, which can be observed in the inelastic neutron 
scattering experiments: 



r 1 = ttt i - 



i 



21n^ 



(23) 



This also agrees with thermal Bethe ansatz calculations 



Let us now compute the nuclear relaxation rates. Nu- 
clear spins are coupled to electron degrees of freedom by 
the magnetic hyperfine hamiltonian: 



H 



HF 



(24) 



I is the nuclear spin, S is the electron spin, a enumer- 
ates spin projections for sites i and j . We will use the 
following expressions p2[ for calculating T\ and T 2 g- 



1 



1 

T20 



2k B T 

h 2 

P 

8h 2 



dq 



Imx(q,oJo) 



CJ 



2n 



dq 
27 



A\{q)x{q) 



(25) 
(26) 



Here Ay {q) and A±(q) are the hyperfine couplings par- 
allel and perpendicular to the easy axis of the crystal, loq 
is the nuclear resonance frequency, which is much smaller 
than any other electron energy scale. The magnetic field 
is directed along the c-axis. The q-dependence is smooth 
and arises from appropriate form factors. The suscepti- 
bility x should, in principle, include contributions from 
both the uniform and staggered spin fluctuations. How- 
ever, simple power counting |p^| shows that the staggered 
component is dominant at small T. For the purpose of 
comparison of our theory with experiment it is conve- 



nient to define normalized dimensionless NMR rates p3[ 
which should be universal functions of T / J: 



hJ 



(VT/T 2 G)norm = 



k B T 
P J 



1/2 



n.i 



(27) 



(28) 



A complete calculation of the NMR relaxation rates 
gives: 



(i/ri), 



T/T 2G ) 



2G /norm 




(29) 
(30) 
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Here D = 1/(2tt) 3 / 2 is the non-universal amplitude, 
7 ~ 0.5772157 is the Euler's constant, while the integrals 
Io and Ii are given by 



h 



h 



POO 

/ dx 

10 


r(i^s) 




/>oo 

/ dx 

10 


r(i^£) 





71.2766 



x i?e 



-259.94, (31) 



where ^(x) is digamma function. The l/ln(A/T) term 
could be incorporated to redefine the cutoff A as in Ref. 
IJ. Thus, up to terms 0(1/ ln 2 (A/T)) the temperature 
dependence for l/T% or \TT jT^Q is actually given by the 
square root of the log and loglog terms in the numerator 
of Eq. (pl|) . The ratio of the relaxation rates, however, is 
only weakly temperature dependent. We find: 



1.680 1 



0.7632 \ 
HA/Tj) 



(32) 



I/T2G, which can be incorporated as a correction to the 
non-universal cutoff scale A. The relaxation rate ratio 
(T2G/T1VT), however, filters out the common factor, and 
therefore picks up weak l/ln(A/T) temperature depen- 
dence, which we have calculated. We note that our result 
is similar to the anzats used by Starykh et al jl3| . There 
are, however, important differences. Starykh et al. |T5| ] 
don't have the loglog term, which turns out to be the 
most important correction in this approximation. The 
1 / ln(To /T) weak temperature dependence for the ratio 
of the relaxation rates was also not explicitly obtained in 
Ref. p3|. Our theoretical results are in excellent agree- 
ment with experimental data of Takigawa et al. 
on Sr 2 Cu03, as shown in Fig.l. 
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FIG. 1. NMR Ti and T 2G /TiT 1/2 vs T/J from Takigawa 



et al. 



fitted to our expression. 



To summarize, the new effect of the higher-order cor- 
rections in the leading irrelevant operator to the dy- 
namic spin susceptibility is, apart from the loglog term 
in the common factor 0n(A/T) + 0.5 ln[ln(A/T)], the 
temperature-dependent anomalous dimension and log- 
arithmic corrections to the correlation length. These 
only lead to an additional weak 0(l/ln(A/T)) tem- 
perature dependence for the relaxation rates \jT\ and 
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